Lack of association between HLA and asymptomatic SARS-CoV-2 infection (preprint)

Human genetic studies of critical COVID-19 pneumonia have revealed the essential role of type I interferon-dependent innate immunity to SARS-CoV-2 infection. Conversely, an association between the HLA-B*15:01 allele and asymptomatic SARS-CoV-2 infection in unvaccinated individuals was recently reported, suggesting a contribution of pre-existing T cell-dependent adaptive immunity. We report a lack of association of classical HLA alleles, including HLA-B*15:01, with pre-omicron asymptomatic SARS-CoV-2 infection in unvaccinated participants in a prospective population-based study in the US (191 asymptomatic vs. 945 symptomatic COVID-19 cases). Moreover, we found no such association in the international COVID Human Genetic Effort cohort (206 asymptomatic vs. 574 mild or moderate COVID-19 cases and 1,625 severe or critical COVID-19 cases). Finally, in the Human Challenge Characterisation study, the three HLA-B*15:01 individuals infected with SARS-CoV-2 developed symptoms. As with other acute primary infections, no classical HLA alleles favoring an asymptomatic course of SARS-CoV-2 infection were identified. These findings suggest that memory T-cell immunity to seasonal coronaviruses does not strongly influence the outcome of SARS-CoV-2 infection in unvaccinated individuals.

Genome-wide Association Study of Long COVID (preprint)

Infections can lead to persistent or long-term symptoms and diseases such as shingles after varicella zoster, cancers after human papillomavirus, or rheumatic fever after streptococcal infections(1,2). Similarly, infection by SARS-CoV-2 can result in Long COVID, a condition characterized by symptoms of fatigue and pulmonary and cognitive dysfunction(3-5). The biological mechanisms that contribute to the development of Long COVID remain to be clarified. We leveraged the COVID-19 Host Genetics Initiative(6,7) to perform a genome-wide association study for Long COVID including up to 6,450 Long COVID cases and 1,093,995 population controls from 24 studies across 16 countries. We identified the first genome-wide significant association for Long COVID at the FOXP4 locus. FOXP4 has been previously associated with COVID-19 severity(6), lung function(8), and cancers(9), suggesting a broader role for lung function in the pathophysiology of Long COVID. While we identify COVID-19 severity as a causal risk factor for Long COVID, the impact of the genetic risk factor located in the FOXP4 locus could not be solely explained by its association to severe COVID-19. Our findings further support the role of pulmonary dysfunction and COVID-19 severity in the development of Long COVID.

Exome-wide association study to identify rare variants influencing COVID-19 outcomes: Results from the Host Genetics Initiative (preprint)

Host genetics is a key determinant of COVID-19 outcomes. Previously, the COVID-19 Host Genetics Initiative genome-wide association study used common variants to identify multiple loci associated with COVID-19 outcomes. However, variants with the largest impact on COVID-19 outcomes are expected to be rare in the population. Hence, studying rare variants may provide additional insights into disease susceptibility and pathogenesis, thereby informing therapeutics development. Here, we combined whole-exome and whole-genome sequencing from 21 cohorts across 12 countries and performed rare variant exome-wide burden analyses for COVID-19 outcomes. In an analysis of 5,048 severe disease cases and 571,009 controls, we observed that carrying a rare deleterious variant in the SARS-CoV-2 sensor toll-like receptor TLR7 (on chromosome X) was associated with a 5.3-fold increase in severe disease (95% CI: 2.75-10.05, p=5.41x10-7). These results further support TLR7 as a genetic determinant of severe disease and suggest that larger studies on rare variants influencing COVID-19 outcomes could provide additional insights.

Evidence for SARS-CoV-2 Delta and Omicron co-infections and recombination (preprint)

Between November 2021 and February 2022, SARS-CoV-2 Delta and Omicron variants co-circulated in the United States, allowing for co-infections and possible recombination events. We sequenced 29,719 positive samples during this period and analyzed the presence and fraction of reads supporting mutations specific to either the Delta or Omicron variant. Our sequencing protocol uses hybridization capture and is thus less subject to artifacts observed in amplicon-based approaches that may lead to spurious signals for recombinants. We identified 20 co-infections, one of which displayed evidence of a low recombinant viral population. We also identified two independent cases of infection by a Delta-Omicron recombinant virus, where 100% of the viral RNA came from one clonal recombinant. In both cases, the 5'-end of the viral genome was from the Delta genome, and the 3'-end from Omicron, though the breakpoints were different. Delta-Omicron recombinant viruses were rare, and there is currently no evidence that the two Delta-Omicron recombinant viruses identified are more transmissible between hosts compared to the circulating Omicron lineages.

Wastewater sequencing uncovers early, cryptic SARS-CoV-2 variant transmission (preprint)

As SARS-CoV-2 becomes an endemic pathogen, detecting emerging variants early is critical for public health interventions. Inferring lineage prevalence by clinical testing is infeasible at scale, especially in areas with limited resources, participation, or testing/sequencing capacity, which can also introduce biases. SARS-CoV-2 RNA concentration in wastewater successfully tracks regional infection dynamics and provides less biased abundance estimates than clinical testing. Tracking virus genomic sequences in wastewater would improve community prevalence estimates and detect emerging variants. However, two factors limit wastewater-based genomic surveillance: low-quality sequence data and inability to estimate relative lineage abundance in mixed samples. Here, we resolve these critical issues to perform a high-resolution, 295-day wastewater and clinical sequencing effort, in the controlled environment of a large university campus and the broader context of the surrounding county. We develop and deploy improved virus concentration protocols and deconvolution software that fully resolve multiple virus strains from wastewater. We detect emerging variants of concern up to 14 days earlier in wastewater samples, and identify multiple instances of virus spread not captured by clinical genomic surveillance. Our study provides a scalable solution for wastewater genomic surveillance that allows early detection of SARS-CoV-2 variants and identification of cryptic transmission.

HLA-A*03:01 is associated with increased risk of fever, chills, and more severe reaction to Pfizer-BioNTech COVID-19 vaccination (preprint)

COVID-19 vaccines are safe and highly effective, but some individuals experience unpleasant reactions to vaccination. As the majority of adults in the US have received a COVID-19 vaccine this year, there is an unprecedented opportunity to study the genetics of reactions to vaccination via surveys of individuals who are already part of genetic research studies. Here, we have queried 17,440 participants in the Helix DNA Discovery Project and Healthy Nevada Project about their reactions to COVID-19 vaccination. Our GWAS identifies an association between severe vaccine side effects and HLA-A*03:01. This association was statistically significant only for those who received the Pfizer-BioNTech vaccine (BNT162b2; p=4.70E-11), but showed a trending association in those who received the Moderna vaccine (mRNA-1273; p=0.005) despite similar sample sizes for study. In Pfizer-BioNTech recipients, HLA-A*03:01 was associated with a two-fold increase in risk of severe vaccine side effects. The effect was consistent across ages, sexes, and whether the person had previously had a COVID-19 infection. The reactions experienced by HLA-A*03:01 carriers were driven by associations with chills, fever, fatigue, and in general feeling unwell.

Rapid displacement of SARS-CoV-2 variant B.1.1.7 by B.1.617.2 and P.1 in the United States (preprint)

The SARS-CoV-2 variant of concern B.1.617.2 displaced B.1.1.7 as the dominant variant in England and other countries. This study aimed to determine whether B.1.617.2 was also displacing B.1.1.7 in the United States. We analyzed PCR testing results and viral sequencing results of samples collected across the United States, and showed that B.1.1.7 was rapidly being displaced and is no longer responsible for the majority of new cases. The percentage of SARS-CoV-2 positive cases that are B.1.1.7 dropped from 70% in April 2021 to 42% in just 6 weeks. Our analysis showed rapid growth of variants B.1.617.2 and P.1 as the primary drivers for this displacement. Currently, the growth rate of B.1.617.2 was higher than P.1 in the US (0.61 vs. 0.22), which is consistent with reports from other countries. Lastly, we showed that B.1.617.2 was growing faster in counties with a lower vaccination rate.

Genomic epidemiology identifies emergence and rapid transmission of SARS-CoV-2 B.1.1.7 in the United States (preprint)

As of January of 2021, the highly transmissible B.1.1.7 variant of SARS-CoV-2, which was first identified in the United Kingdom (U.K.), has gained a strong foothold across the world. Because of the sudden and rapid rise of B.1.1.7, we investigated the prevalence and growth dynamics of this variant in the United States (U.S.), tracking it back to its early emergence and onward local transmission. We found that the RT-qPCR testing anomaly of S gene target failure (SGTF), first observed in the U.K., was a reliable proxy for B.1.1.7 detection. We sequenced 212 B.1.1.7 SARS-CoV-2 genomes collected from testing facilities in the U.S. from December 2020 to January 2021. We found that while the fraction of B.1.1.7 among SGTF samples varied by state, detection of the variant increased at a logistic rate similar to those observed elsewhere, with a doubling rate of a little over a week and an increased transmission rate of 35-45%. By performing time-aware Bayesian phylodynamic analyses, we revealed several independent introductions of B.1.1.7 into the U.S. as early as late November 2020, with onward community transmission enabling the variant to spread to at least 30 states as of January 2021. Our study shows that the U.S. is on a similar trajectory as other countries where B.1.1.7 rapidly became the dominant SARS-CoV-2 variant, requiring immediate and decisive action to minimize COVID-19 morbidity and mortality.

S gene dropout patterns in SARS-CoV-2 tests suggest spread of the H69del/V70del mutation in the US. (preprint)

Recently, multiple novel strains of SARS-CoV-2 have been found to share the same deletion of amino acids H69 and V70 in the virus S gene. This includes strain B.1.1.7 / SARS-CoV-2 VUI 202012/01, which has been found to be more infectious than other strains of SARS-CoV-2, and its increasing presence has resulted in new lockdowns in and travel restrictions leaving the UK. Here, we analyze 2 million RT-PCR SARS-CoV-2 tests performed at Helix to identify the rate of S gene dropout, which has been recently shown to occur in tests from individuals infected with strains of SARS-CoV-2 that carry the H69del/V70del mutation. We observe a rise in S gene dropout in the US starting in early October, with 0.25% of our daily SARS-CoV-2-positive tests exhibiting this pattern during the first week. The rate of positive samples with S gene dropout has grown slowly over time, with last week exhibiting the highest level yet, at 0.5%. Focusing on the 14 states for which we have sufficient sample size to assess the frequency of this rare event (n>1000 SARS-CoV-2-positive samples), we see a recent expansion in the Eastern part of the US, concentrated in MA, OH, and FL. However, we cannot say from these data whether the S gene dropout samples we observe here represent the B.1.1.7. strain. Only with an expansion of genomic surveillance sequencing in the US will we know for certain the prevalence of the B.1.1.7 strain in the US.

Using symptom-based case predictions to identify host genetic factors that contribute to COVID-19 susceptibility (preprint)

Epidemiological and genetic studies on COVID-19 are hindered by inconsistent and limited testing policies to confirm SARS-CoV-2 infection. Recently, it was shown that it is possible to predict potential COVID-19 cases using cross-sectional self-reported disease-related symptoms. Using a previously reported COVID-19 prediction model, we show that it is possible to conduct a GWAS on predicted COVID-19 which benefits from a larger sample size in order to gain new insights into the genetic susceptibility of the disease. Furthermore, we find suggestive evidence that genetic variants for other viral infectious diseases do not overlap with COVID-19 susceptibility and that severity of COVID-19 may have a different genetic architecture compared to COVID-19 susceptibility. Our findings demonstrate the added value of using self-reported symptom assessments to quickly monitor novel endemic viral outbreaks in a scenario of limited testing. Should there be another outbreak of a novel infectious disease, then we recommend repeatedly collecting data of disease-related symptoms.